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ABSTRACT 

We report on 12 years of observations of PSR J1713+0747, a pulsar in a 68-day orbit with a white 
dwarf. Pulse times of arrival were measured with uncertainties as small as 200 ns. The timing data 
yielded measurements of the relativistic Shapiro delay, perturbations of pulsar orbital elements due 
to secular and annual motion of the Earth, and the pulsar's parallax, as well as pulse spin-down, 
astrometric, and Keplerian measurements. The observations constrain the masses of the pulsar and 
secondary star to be mi = 1.3 ± 0.2 M Q and m 2 = 0.28 ± 0.03 M , respectively (68% confidence). 
Combining the theoretical orbital period-core mass relation with the observational constraints yields a 
somewhat higher pulsar mass, mi = 1.53lQQg M©. The parallax is n = 0.89±0.08mas, corresponding 
to a distance of 1.1 ± 0.1 kpc; the precision of the parallax measurement is limited by uncertainties in 
the electron content of the solar wind. The transverse velocity is unusually small, 33 ± 3kms -1 . We 
find significant timing noise on time scales of several years, but no more than expected by extrapolating 
timing noise statistics from the slow pulsar population. With the orientation of the binary orbit fully 
measured, we are able to improve on previous tests of equivalence principle violations. 
Subject headings: stars: neutron — binaries: general — pulsars: individual (PSR J1713+0747) 



1. INTRODUCTION 

High precision timing of a radio pulsar binary system 
can reveal a wealth of information about the dynamics of 
the binary, the astrometry of the system, and the natures 
of the pulsar and secondary stars. PSR J1713+0747, a 
4.6 ms pulsar in a 67.8 day orbit with a white dwarf, is 
among the very best pulsars for high precision timing: 
it has a hig h flux density, sha llow spectrum, and sharp 
pulse peak ({Foster et alJH993|) . 

We observed PSR J1713+0747 over six years as part 
of a systematic high precision pulsar timing study at the 
Arecibo Observatory. We used data acquisition systems 
employing coherent dedispersion, allowing substantial 
improvements in timing precision compared to previous 
work. 

Early obse r vation s of this source were reported by 
iCamilo et"al1 l|1994fl (herein CFW), who analyzed 1.5 
years of pulsar timing data. Previous observations have 
also been reported bv lvan Straten fc Bailesl l)2003lh We 
combined the data of CFW with our newer observations 
to produce a single data set spanning twelve years. 
The superior timing precision of the new data and the 
longer total time span of observations yield substantial 
refinements of all measurements reported by CFW, as 
well as a number of additional new measurements. 



High precision detection of the Shapiro delay allows 
the pulsar and secondary star masses to be separately 
measured. The Shapiro delay measurement reported in 
CFW did not have sufficient precision to independently 
determine the pulsar and secondary star masses. These 
quantities are of interest because, while measured double 
neutron s tar masses fall into a very n arrow range, 1.25 to 
1.44 M© ijLattimer fc Prakashll2004[K extended accretion 
during the formation of pulsar-white dwarf binaries may 
lead to higher pulsar masses in these systems. 

High precision timing of pulsars can be used to 
place limits on the gravitational wave background at 
frequencies of 10~ 6 to 10 Hz by searc hing for arrival 
time variations on tim e scales of years l|Lommen et alJ 
120031 iKaspi et al] |l994j), but such measurements require 
that the timing signal not be contaminated by "timing 
noise," random variations of the pulsar rotation. 
Timing nois e is inversely correla t ed wi th rotation period 
derivative ijArzoumanian et al.l Il994j) . so millisecond 
pulsars, which have very small period derivatives, are 
prime candidates for gravitational wave studies. The 
12-year time span of data on PSR J1713+0747 has 
revealed significant pulse arrival variations beyond those 
expected from simple magnetic dipole spin-down. We 
quantify this apparent timing noise, but we leave the 
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analysis of the PSR J1713+0747 signal in the context of 
the gravitational wave background to another work. 

We detected both annual and secular perturbations 
of orbital elements due to the changing Earth-binary 
line-of-sight. This allowed us to uniquely determine 
the inclination and position angle of the orbit. This 
is only the second pulsar binary system for which 
timing observations have full y specified both angles 
of the orbital orientation (van Straten et all I200H 
Ivan Straten fc Bailesll200^ . 

Preliminary results from this project, using data 
collec ted throug h mid- 2002, were publ ished elsewhere 
fe.g.. iNice et"aLll200l INice et all 12004]) . The present 
paper incorporates two further years of high quality data, 
resulting in some changes to the best-fit parameters, but 
within the expected uncertainties. 

The plan of this paper is as follows. In §2 we summarize 
the observations. In §3 we describe the timing model. In 
§4 we present the parameters of the pulsar and binary 
system derived from fitting the timing model to the 
observational data. In §5 we analyze the timing model 
parameters in the contexts of timing noise, stellar masses, 
pulsar distance and velocity, and theories of strong field 
gravitation. In §6 we summarize the key results. 

2. OBSERVATIONS 

The 305-m Arecibo radio telescope recorded a total of 
343 pulse times of arrival (TOAs) of PSR J1713+0747 
on 166 separate days between April 1992 and May 2004. 
Table ^ summarizes the observations. A shutdown of 
the telescope for an upgrade resulted in a gap in the 
observations between 1994 and 1998. 

Observations made in 1992 and 1993 are described in 
CFW. Briefly, those data consist of 68 TOAs collected at 
intervals of about t wo weeks using the Prin ceton Mark III 
observing system (Stincbri ng et alJll992[) coupled to a 
2 x 32 x 1.25 MHz filter bank, with 32 spectral channels for 
each sense of polarization. In every scan, the square-law 
detected outputs of each pair of channels (representing 
two polarizations) were smoothed with time constants 
of 20 fj,s or 78 /us, summed, folded synchronously at the 
predicted pulse period, and shifted in time to remove 
dispersive delay. The channels were then summed to 
create total-intensity pulse profiles of 256 bins. 

From 1998 through 2004, ob servations were m ade 
using Princeton Mark IV system ijStairs et al.ll2000j) and 
the Arecibo-Berkeley Pulsar Processor (ABPP), often 
running in parallel to analyze the same radio frequency 
signal. A total of 275 TOAs were obtained on 132 
different days, usually at intervals of a few weeks, 
but occasionally more densely sampled. A typical day 
included an hour of observations at 1410 MHz and a half 
hour of observations at 2380 MHz. 

The Mark IV system critically samples and records 10 
MHz pass bands in each sense of circular polarization, 
quantized with 2-bit resolution. The recorded voltages 
are analyzed off line: the data stream is coherently 
dedispersed, after which the self- and cross-products of 
the voltages are calculated and folded synchronously 
at the pulse period. Observations are continuous over 
(typically) 29 minute intervals, but are analyzed in blocks 
of 190 seconds, each of which yields a 1024 bin pulse 
profile with four Stokes parameters. 

The ABPP filters the passband into narrow spectral 



channels, samples voltages with 2-bit resolution, and 
applies coherent dedispersion to each channel using 
3-bit coefficients. For PSR J1713+0747 at 1410 MHz, 
thirty- two spectral channels of width 1.75 MHz are 
processed in each polarization, for a total bandwidth 
of 56 MHz; at 2380 MHz, the channel bandwidths are 
3.5 MHz, for a total bandwidth of 112 MHz. The 
dedispersed time series in each channel is folded 
synchronously at the pulse period and integrated for 180 
seconds. 

We used conventional techniques to measure pulse 
arrival times. For the Mark III and Mark IV data, 
each pulse profile was cross-correlated with a standard 
template to measure the phase offset of the pulse within 
the profile. For the ABPP data, an analytical model 
of the pulse profile was fit to each data profile. In 
either case, the time offset measured from the profile was 
added to the start time of the integration and translated 
to its middle to yield a TOA. The start times were 
referenced to the observatory atomic clock, which was 
corrected retroactively to the UTC timescale using data 
from the Global Positioning System (GPS) satellites. 
For each data acquisition system, all TOAs at a given 
frequency collected on a given day were averaged to make 
a single effective TOA. The time intervals spanned by 
these average points is listed as "typical integration" in 
Table [B 

The ABPP and Mark IV systems often ran in parallel, 
analyzing the same radio frequency signal processed 
through the same amplifiers and many of the same 
filters. Of the TOAs summarized in Tabled seventy-six 
pairs of ABPP and Mark IV TOAs were collected 
simultaneously. Nevertheless, we have treated them as 
independent data streams in the timing analysis. There 
are two justifications for this approach. First, the width 
of the passband measured by the ABPP is substantially 
larger than that measured by Mark IV, so the signal 
measured by the two machines are somewhat different. 
Second, empirical tests of the residual arrival times, after 
removing the best fit timing model, show only modest 
correlation between TOAs measured by the two systems. 
Correlation coefficients are 0.33 and 0.30 for the 1410 and 
2380 MHz residual TOAs, respectively. 

We found that formal measurement uncertainties 
calculated directly from the TOA measurements 
of individual data records tended to moderately 
underestimate the true scatter in the arrival times. The 
cause of the underestimation is not known, but it is 
a common phenomenon in millisecond pulsar timing. 
A particular challenge for these observations is the 
correction of the coarse quantization of incoming signals 
given the highly variable nature of the signals. In any 
case, we added "systematic" terms in quadrature to these 
formal uncertainties in order to produce timing fits with 
reduced \ 2 values close to 1. The root-mean-square 
(RMS) residual arrival times listed in Table ^ were 
calculated by giving all TOAs equal weight, independent 
of uncertainties. 

3. TIMING MODEL 

We used the tempo 1 software package to fit a 
pulse timing model to the observed TOAs. The 

1 http:/ /pulsar. princeton.edu/tcmpo 
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TABLE 1 

Summary of Observations 



System Dates Frequency Bandwidth Number Typical RMS 

(MHz) (MHz) of Integration Residual" 

TOAs (min) Qjb) 



Mark III 1992.6-1993.0 1400 

1992.3-1994.1 1400 

ABPP 1998.1-2004.4 1410 

1999.7- 2004.4 2380 
Mark IV 1998.6-2004.4 1410 

1999.8- 2004.4 2380 



40 9 47 0.70 b 

40 59 47 0.41 c 

56 101 60 0.18 

112 49 30 0.35 

10 81 58 0.28 

10 44 29 0.46 



a Values incorporate the effect of averaging TOAs from shorter integration times. 
b Filter bank used a 78-fis time constant. 
c Filter bank used a 20-fis time constant. 



model incorporated pulsar rotation, astrometry, orbital 
motion, and dispersion of the pulsar signal by the 
ionized interstellar medium. The fits allowed for 
arbitrary time offsets between sets of TOAs taken 
with different observing systems and at different 
frequencies. Earth m otion was modeled u sing the 
JPL DE405 ephemeris ijStandisbl Il99l l200l see also 
|03>- The ultimate time reference was the TT(BIPM03) 
scale {Bureau I nternati onal des Poids et Mesuresl 2003) 
adjusted to bar ycenter time scale TDB w ith the TE405 
time ephemeris ijlrwin fc Fukus hima 1992). 

In this section, we describe the elements of the timing 
model in detail. We defer discussion of the results of the 
model fit until 21 However, we will make reference of 
these results to justify components of the timing model. 




1995 2000 
Year 



Fig. 1— Post-fit residuals of PSR J1713+0747. (a) Pulsar 
rotation modeled by vq and v\ only; timing noise is evident. 
Because the early data (pre-1995) are less precise than the later 
data, they are down-weighted in the fit; as a result, the fit allows 
their residual arrival times show greater variation. Alternate values 
of uq and v\ could be found which reduced the timing noise in the 
early data at the expense of increasing it in the later data, (b) 
Pulsar rotation modeled by vq through ug. Error bars shown here 
do not include compensation for systematic uncertainties, although 
such compensation was included in these timing models. 



3.1. Pulsar Rotation and Timing Noise 

The pulsar rotation frequency at time t can be 
written as a polynomial expansion, i>{t) = Vq + 
v\(t — to) + \ v 2{t — ^o) 2 + where t is a reference 
epoch near the center of the data span and vq, V\, 
V2, etc., are the pulsar rotation frequency, frequency 
derivative, frequency second derivative, and so on. 
Under the standard model of a neutron star with a 
rotating magnetic dipole, a millisecond pulsar such as 
PSR J1713+0747 should have negligibly small values 
of V2 and higher order terms. However, timing fits 
incorporating only vq and V\ leave a systematic signature 
in the residual arrival times after removing the best-fit 
model (Figure [T]l. To remove the long-term systematic 
trend, we "whitened" the data by incorporating seven 
additional terms to the model, i>2 through i>g. These 
terms correspond to time scales on the order of years, 
much longer than the orbital period of the pulsar, and 
hence they have negligible effect on the measurement of 
orbital elements. 

3.2. Dispersion Measure Variations 

Radio pulses traversing the solar system and the 
interstellar medium are delayed by dispersion. The time 
delay in seconds is At BM = DM /(2.41 x 10~ 16 / 2 ), 
where / is the radio frequency in Hz and the integral of 
electron density along the line of sight, DM = J n e (l) dl, 
is the dispersion measure in peern - . Observations at 
two frequencies, 1410 and 2380 MHz, over the last 5 
years of observations, allow us to search for variations 
in DM over time. Of particular concern are variations 
due to propagation through the ionized solar wind as 
the Earth moves about the Sun. We modeled the 
solar electron density at distance r from the Sun as 
n e (r) = rto(lAu/r) 2 , where no is the electron density 
at r=lAU which we take as a free parameter in the 
pulsar timing solution. The solar contribution to DM 
is calculated by integrating n e (r) along the path from 
Earth toward the pulsar. Ulysses data have shown the 
1/r 2 scali ng to hold over a wi de range of heliocentric 
latitudes <|Issautier et all 12001). The 1/r 2 model is 
clearly an oversimplification — it neglects the substantial 
difference between the high density slow wind along the 
ecliptic and the low density wind at high latitudes, and 
it does not allow for temporal variations in no. However, 
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because it has a simple analytic form and requires only a 
single additional parameter in the fit, it is a convenient 
form to use for the timing model. The best fit electron 
density at 1 AU is no = 5 ± 4 electrons per cm 3 . This 
gives peak-to-peak DM variations of 0.0002 pc cm -3 , and 
arrival time variations of 400 ns at 1410 MHz. 

3.3. Orbital Kinematics Including Shapiro Delay 

Orbital kinematics were incorporated into the 
timing model by means of the theory-independent 
representation of Damour & Deruelle (1986). Parameters 
of the orbital model include (1) five Keplerian orbital 
elements: orbital period, semi-major axis projected 
into the line of sight, x = (ai sini) / c, where a\ is the 
semi-major axis, i the inclination angle, and c the speed 
of light; eccentricity, e; angle of periastron, u>; and time 
of periastron passage, T ; (2) secular variations of the 
Keplerian elements, most notably the time derivative of 
x, denoted x; (3) the orientation of the system, defined 
by the inclination of the orbit, i, where 0° < i < 180°, 
and the position angle of ascending node, f2, where 
0° < ft < 360° and ft is defined north through east 2 
(see Figure EJ); and (4) the masses of the pulsar and the 
secondary star, mi and m-i- 

The masses and inclination are connected by the mass 
function. 



_ (m 2 sini) 3 
Ji = 



(1) 



x 6 (2ir 

(mi + m 2 ) 2 T \P b/ 

where T = GM /c 3 = 4.925 x 10 _6 s. In the analysis 
below, we treat i and mi as independent parameters in 
the timing model, and we use equation ^ to determine 
m\. 

According to general relativity, pulses are retarded as 
they propagate through the gravitational potential well 
of the secondary. For a nearly circular orbit, this Shapiro 
delay is 



At s = — 2 TO2 Tq ln[l — sin i sin( 



■ w ) 



(2) 




Tangent Plane of the Sky 



Fig. 2. — Geometry of the binary system, showing position 
angle of ascending node, f2, inclination, i, and angle of periastron, 
Lo. The black dot indicates the ascending node. 

2 For celestial position angles, we follow the standard convention 
that 0° is north and 90° is east. For inclination, we follow the 
convention that an orbit with i = 0° has an angular momentum 
vector pointing toward E arth. These definitions differs from 
Kopeikin (1995, 1996) and Ivan Straten et ahl 12001]) . who define 
Q such that 0° is east and 90° is north, and define i = 0° to be an 
orbit with angular momentum vector pointing away from Earth. 



where <f> is the orbital phase. In principle, measurement 
of the Shapiro delay yields mi and sini. In practice, 
unless sin i ~ 1, the Shapiro delay is highly covariant 
with x in the timing fit and hence difficult to measure. 
Nevertheless, as shown in Figure the Shapiro delay 
can be clearly distinguished in the PSR J1713+0747 
data. The constraints on 7712 and i that arise from this 
measurement are quantified in 21 Since 0° < i < 180° 
(see Figure |2J), there in an ambiguity in deriving i from 
sini in equation [21 because both i and 180° — i are 
solutions. The resolution of this ambiguity is discussed 
below. 

3.4. Projection Effects due to Proper Motion 

Proper motion o f a binary syst em results in secular 
changes in u> and x <|Kopeikinll996j) . For a nearly circular 
orbit, a secular change in to is indistinguishable from a 
small perturbation of the orbital period, and hence is 
unmeasurable. In contrast, the secular change in x is 
significant. The time derivative ol x, x, is given by 



H cot i sin(0 M — ft), 



(3) 



where \i and # M are the magnitude and position angle of 
the proper motion, respectively. In effect, i is determined 
by the Shapiro delay f M3.3[l . and ft is then constrained 
by the x measurement. Solving equation |21 for f2 yields 
two possible values for each of the two values of i 
f H3.3l) . Thus there are four distinct binary orientations 
(combinations of i and f2) allowed by the Shapiro delay 
and x measurements. 

3.5. Annual- Orbital Parallax 

The Earth's annual motion about the solar system 
barycenter changes the line-of-sight to the binary system 




1 ]_ (c) J J I j j J T j 
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Fig. 3. — Shapiro delay in the timing residuals of 
PSR J1713+0747 as a function of orbital phase, (a) The Shapiro 
delay term is excluded from the model while the rest of the 
parameters are fixed at their best-fit values from Tablel2l The solid 
curve shows the general relativistic delay predicted by equation |3 
(b) The Shapiro term is excluded from the model but the other 
parameters are permitted to vary, (c) Shapiro delay is included in 
the model. The residuals in this panel are from the same timing 
fit as those in Figure lib. 



5 



and perturbs the binary parameters x and to, an 
effect known as annual- orbital parallax (|KoDeikiiJ|1995t 
Ivan Straten et alJl2001t Ivan Straten fc Bairesll2003jh We 
incorporated annual-orbital parallax into the timing 
model by perturbing the values of x and uj before 
calculating pulse arrival time delays using a standard 
orbital model. The perturbation formulae are given by 
Kopeikin (1995); we summarize them here. The position 
of the Earth relative to the solar system barycenter is 
given by the vector X, Y, Z. The pulsar position is right 
ascension a and declination S. Define 



VJO 



Aio = —X sin a + Y cos a, 
—X sin S cos a — Y sin S sin a 



Z cos 5. 



(4) 



These are the east and north motions of the Earth in a 
coordinate system parallel to the plane of the sky. The 
observed values of x and uj are perturbed from their 
intrinsic values by 



•^obs — *£int 



and 



1 



-"int 



cot i 



(Aio cos Vt — Ajq sin 0) 



(5) 



CSC I 



(Aio sin f2 + A jo cos 17). (6) 



We found the annual-orbital perturbation of x to have 
only a marginal effect on the timing of PSR J1713+0747, 
and no useful measurements can be derived from it. 
(Nevertheless, the x perturbation was included in the 
timing model, since it introduced no additional degrees of 
freedom.) On the other hand, the perturbations of uj are 
sufficiently large that incorporating them into the timing 
model significantly improves the goodness of the timing 
fit. FigureEJshows the x 2 values for fits with and without 
the annual-orbital parallax perturbations. In these fits, 
mi, i, x, and f2 were treated as independent parameters, 
unconstrained by equation El hence the best-fit range of 
f2 in Figure 0] differs somewhat from that given in §4. 
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Fig. 4. — The effect of annual-orbital parallax on the timing 
solution is shown by displaying the difference in goodness of fit, 
Ax 2 , as a function of position angle of ascending node. The 
value A\ 2 = corresponds to the best timing model without 
annual-orbital perturbations. The 95% confidence range of f2 is 
indicated. 



As a practical matter, for PSR J1713+0747 the 
annual-orbital parallax does little to improve the 
precision of the timing parameter measurements (since 
i and SI are better measured by Shapiro delay and x). 
However, it allows the fourfold ambiguity in i and £1 to 
be broken, picking out one distinct orientation of the 
orbit. 

4. TIMING ANALYSIS: PARAMETER VALUES 

We fit the measured TOAs to a timing model 
incorporating the phenomena described in fEl We used 
a hybrid procedure to determine the timing parameters 
and their uncertainties. Standard least-squares methods 
are adequate for fitting most of the quantities in the 
timing model. However, for the orientation and mass 
parameters, i, f2, and m%, the \ 2 surfaces are not 
ellipsoidal in the parameter space of interest. To 
investigate the allowed ranges of these parameters 
(and, ultimately, the other parameters as well), we 
analyzed timing solutions over a uniformly sampled three 
dimensional grid of trial values of cosi, m%, and Q, in 
the vicinity of the x 2 minimum. For each combination 
of cosz, mi-, an d ^, we calculated x, the Shapiro 
delay parameters, and the annual-orbital perturbation 
corrections, according to equations [2 El El an d El We 
then performed a timing fit, in which these quantities 
held fixed while all other parameters were allowed to 
vary. We recorded the resulting values of x 2 for each 
grid point. 

The minimum \ 2 timing solution is at cosi = 0.31, 
m 2 = O.28M , and Q = 87°. The parameters 
corresponding to this solution are summarized in TableEl 
Our resu lts show good agreement wi th the less precise 
results of Ivan Straten fc Bailesl 1)20031) . 

We used a statistically rigorous procedure to calculate 
probability distribution functions (PDFs) for the 
individual quantities cosi, m2, and mi and their 
corresponding uncertainties in table El The procedure 
is a straightforward three-dimensional extension of 
the Bayesian algorithm described in appendix A of 
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Fig. 5. — Probability distribution functions for values of fi, mi, 
m2, and cosi. 
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TABLE 2 
Timing Model Parameters' 1 



Measured Quantities 



Right Ascension, a (J2000) 17 h 13 m 49f 5305335(6) 

Declination, <5 (J2000) +07°47'37'.'52636(2) 

Total proper motion, /i (mas yr -1 ) 6.297(7) 

Position angle of proper motion, 6 M 128?66(7) 

Parallax, n (mas) 0.89(8) 

Rotation frequency, v Q (s" 1 ) 218.8118439157321(3) 

First derivative of u , Vi (s -2 ) -4.0835(2) X 10 -16 

Epoch, t (MJD [TDB]) 52000.0 

Dispersion Measure, DMo (pccm -3 ) 15.9960 

Orbital period, P b (days) 67.8251298718(5) b 

Projected semi-major axis, x (lt-s) 32.34242099(2) b 

Eccentricity, e 0.0000749406(13) b 

Time of periastron passage, T (MJD [TDB]). . . 51997.5784(2) b 

Angle of periastron, ui (deg) 176.1915(10) b 

Cosine of inclination angle, cosi 0.31(3) 

Position angle of ascending node, Q (deg) 87(6) 

Companion mass, m2 (Mq) 0.28(3) 



Measured Upper Limits 

First derivative of DM, DMi (pc cm-Syr -1 ) . . . 0(2) X 10~ 5 

First derivative of u), Co (degyr -1 ) 6(4) X 10~ 4 

First derivative of P b , P b 0(6) X 10 _1: 



Derived Quantities 

Proper motion in a, fi a = a cos 6 (mas yr -1 ) . . . 4.917(4) 

Proper motion in <5, fig = 5 (mas yr — 1 ) —3.933(10) 

Galactic longitude, I 28? 748 

Galactic latitude, b 25? 222 

Distance, d (kpc) 1-1(1) 

Rotation period, P (s) 0.004570136525082781(6) 

Observed rotation period derivative, P bs 8.5288(3) X 10 -21 

Intrinsic rotation period derivative, Pi nt 8.1 X 10 -21 

Characteristic age, r (yr) 8 X 10 9 

Magnetic field, B (G) 1.9 x 10 s 

Mass function, /i (Mq) 0.0078962167 

Pulsar mass, m\ (M©) 1-3(2) 

First derivative of x, x (10 -15 ) 6.7(2) 



a Figures in parentheses are 68% confidence uncertainties in the last digit 
quoted. 

b Kcplcrian orbital elements P b , x, e, To, and lo are covariant with 
post-Keplerian elements cosi, Q, and ni2- This covariance is not reflected in 
the values of Keplerian elements and uncertainties quoted in this table, which 
were derived in a timing model with cosi, f2, and rri2 fixed to their best-fit 
values. 



ISplaver et all lEffll . We assigned a probability to each 
grid point based on the difference between its \ 2 an d the 
minimum % 2 on the grid. After suitable normalization, 
we summed the probabilities of all points associated with 
a given range of mi (or 1112, cosi, or tt) to calculate the 
PDF. 

In effect, we are incorporating a uniform prior 
distribution in cosi, Q, and m,2- The uniform 
distributions in cos i and Q correspond to the probability 
distributions for the angular momentum vectors of 
randomly oriented orbits, while the uniform distribution 
of mi is simply an ad hoc assumption. 

The calculated PDFs are given in Figure They 
yield the 68% confidence estimates given in Table|2 most 
notably m x = 1.3 ± 0.2 M Q and m 2 = 0.28 ± 0.03 M . 
The 95% confidence estimates on the masses are m\ = 
1.3±^M Q and m 2 = 0.28t^M Q . 



The interdependence between these parameters is show 
in the multidimensional confidence plots in Figures 
and To produce Figure H3 we summed normalized 
probabilities across all values of £1 to yield a two 
dimensional probability distribution in cosi and mi- 
The contours enclose 68% and 95% of the probability 
distribution in this space. In Figure [7| we show 
representative "slices" in constant il of the regions 
enclosing 68% and 95% of the three dimensional grid. 

5. DISCUSSION 

5.1. Orbital Period-Core Mass Relation 

5.1.1. Testing the Relation 

Binary evolution theory predicts a specific relationship 
between the orbital period, Pb, of an evolved neutron 
star-white dwarf binary and the secondary mass, mi, 
which is presumed to be equal to the helium core of 
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the progenitor of the secondary. Measured values of 7712 
can be used to test this relation. Figure |H1 shows the 
rel ations derived from b inary evolut ion tracks calculated 
bv ITauris fe Savoniid (|1909f) and iPodsiadlowski et alJ 
(2002). The la tter curves are the me dian and upper 
values of ni2 of Rappa port et al.1 lll995l). wh i ch pro vide 
good bounds on the tracks of iPfahl et all <)2002|) (P. 
Podsiadlowski, private communication). The models are 
in good agreement with the measured values. 

5.1.2. Implications of the Relation for PSR ,11713+0747 
System Masses 

If, rather than using timing measurements to test 
the orbital period-core mass relation, we accept this 
relation as correct, than we can use it to refine the range 
of pulsar mass allowed by t he timing data. For the 
orbital period of ,71713+0747. ITauris fc Savorniel {,1999) 
calculate the secondary mass to be 0.31 M© < 777,2 < 
0.34 M Q , with the lower end of this range for a population 
I progenitor and the hi g her e nd for a population II 
progenitor. IPfahl et alJ ll200l give a similar range, 
0.30 M Q < to 2 < 0.35 Mq (again usi ng the higher end 
of the range of lRappaport et aLlll995j) . 

Repeating the statistical analysis of timing solutions 
on a cos 7 — m,2 — £1 as in § but restricting 777,2 to 
the range 0.30 M Q < 7772 < 0.35 M Q yields a value of 
mi = l-53io;°6 (68% confidence). 

5.2. Neutron Star Mass 

Masses of pulsars and secondary stars in neutron 
star-neutron star binaries all fall within the range 1.25 
to 1.44 Mq. This is near the minimum mass for neutron 
star formation; the maximu m mass is not known, and 
may range up to 3 M Q (see lLattimer fc Prakashl 12004 
for a review). Because binary systems which evolve into 
pulsar-white dwarf binaries like the PSR J1713+0747 
system undergo extended periods of mass transfer, the 
pulsars in these systems might be expected to be heavier 
than those in neutron star neutron star binaries. 
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iThorsett fc Chakrabartvl l)1999f) studied the entire 
binary pulsar population and found it to be consistent 
with a narrow Gaussian distribution with mean 1.35 Mq 
and width 0.04 M©. Their calculation included a 
statistical analysis of the pulsar-white dwarf population, 
under the assumption that the binaries were randomly 
oriented in space and that the orbital period-core mass 
relation held, and they found that these pulsars could 
be drawn from the same mass distribution as pulsars in 
neutron star-neutron star binaries. 

The measured mass of PSR J1713+0747 based on 
observations alone, mi = 1.3 ± 0.2 M©, is in excellent 
agreement with the 1.35 M© value. In contrast, the mass 
derived when the orbital period-core mass constraint 
is imposed, 7771 = 1.531q'q 6 M©, implies a significantly 
heavier neutron star, perhaps having accreted ~ 0.2 M© 
if initially formed at 1.35 M©. There are two other 
pulsar-helium white dwarf binaries with well measured 
Shapiro delays, PSRs J0437-4715 and B1855+09, which 
have pulsar masses 1.58 ± 0.18M© and L5O±gJ|M , 
respectively ijvan Straten et alJ200HlKaspi et al.ll99|) 3 . 
On the other hand, PSR J2019+2425, in a similar 
system, shows a lack of a detectable Shapiro delay and 
observed secular changes which imply an upper limit 
of 1.51 M,D and a median likelihood value of 1.33 M© 
l)Nice et al.l 12001^1 . The uncertainties on all of these 
measurements are frustratingly large; the question of the 
distribution of pulsars masses in these systems remains 
open. 

5.3. Parallax, Distance, and Velocity 

Because the solar wind introduces annual 
perturbations on the pulse arrival times, the solar 
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Fig. 7. — Contours of 68% and 95% confidence intervals in 
cos i — m'2 — Q space. Four slices of constant Q are shown; in 
each, the intersection of the slice with the 68% and 95% confidence 
region is given. The slice at Q = 95° intersects only the 95% 
region. The dashed line corresponds to the 95% confidence region 
after marginalizing over all values of Q (cf. Figure 0. 



Fig. 6. — Contours of 68% and 95% confidence intervals in 
cos i — rri2 space. 



3 A separate analysis of some of t he Kaspi et al. (1994) 
data by Thorsctt & Chakrabartv 1 1999) reported a smaller mass, 
1.40 ± 0.10 M Q . 



wind electron density parameter, no, is highly covariant 
with parallax, n, and position, a and 6, in the pulse 
arrival time model. For tt, the measurement uncertainty 
using any fixed solar wind model is ±0.05 mas, while 
the uncertainty due to the poor constraint on no 
is ±0.06 mas. We combine these uncertainties in 
quadrature to find tt — 0.89 ± 0.08 mas. As a check, 
a special timing fit was done incorporating only days 
on which multifrequency observations were made, and 
fitting for a separate value of DM on each day on 
which observations were made, yields a parallax of 
0.93 ± 0.24, consistent with our preferred value. The 
parallax measurement corresponds to a distance of 
d= l.l±0.1kpc. 

For PSR J1713+0747's measured value of DM, 
the NE2001 model o f the Galactic electron density 
llCordes fc Laziol 12002^ gives a distance of 0.9 kpc, in 
good agreement with the parallax measurement. 

Right ascension and declination are also highly 
covariant with no; their uncertainties quoted in Table El 
were calculated in the same manner as the parallax 
uncertainty. 

It is well established that millisecond pulsars have 
substantially smaller velocities than the bulk pulsar 
populatio n, with a mean transverse velocity of 85 ± 
lSkms" 1 llCordes fc Chernofi1ll997tlToscano et"allll999t 
Nice fc Tavlorl 11995ft . The transverse velocity of 
PSR J1713±0747, V = fid = 33 ± Skms^ 1 , is small 
even by the standards of millisecond pulsars. 

5.4. Kinematic Corrections to Spin Period Derivative; 
Age & Magnetic Field 

The observed pulse period derivative, P bs = —vi/vq, 
is biased away from its intrinsic v a lue, P int, by Doppler 
accelerations. iDamour fc Tavlorl l)1991j) analyzed this 
bias for orbital period derivatives; their work applies 
equally well to spin period derivatives. The observed 
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Fig. 8. — Theoretical rri2 — Pf, relation and measured white 
dwarf masse s in pulsar— white d warf binaries. Measured mas ses are 
from Ivan Straten et alj l l20ulT) (J0437-4715), IKaspi et all Il99l) 
(B1855+09), and this work (J1713+0747). 



and intrinsic quantities are related by 

-Pint = -fobs — AP PM — AP rot — AP Z , (7) 

where APpm = P[J?d/c is the bias due to proper motion; 

AP ro t = P(&rot • n )/c is the bias due to the relative 
acceleration of pulsar and Earth in the Galactic plane, 
&rot, projected into a unit vector pointing from Earth to 
pulsar, n; and AP Z = P(a z • n)/c is the bias due to the 
acceleration of the pulsar toward the galactic disk, a z , 
projected into the line-of-sight. 

For PSR J1713±0747, we find AP PM = 0.46 x 10" 21 , 
AP rot = 0.22 xl0~ 21 and AP Z = -0. 32xlQ- 21 , the latter 
calculated using the potential of iKuiiken &: Gilmorel 
(1989). The net bias is AP = 0.44 x 10" 21 , so that 
we estimate the intrinsic period derivative to be Pj n t = 
8.1 x 10~ 21 . 

The characteristic age of PSR J1713+0747 is t = 
P/2P; nt = 8 Gyr. This probably overestimates the 
true age of the pulsar. Hansen & Phinney (1998) 
obtain a considerably lower estimate by using the optical 
measurements of Lundgren, Foster, & Camilo (1996) to 
infer the age of the white dwarf to be 6.3 Gyr < r < 
6.8 Gyr. 

A discrepancy between the characteristic age and the 
true age of a millise cond pulsar is not uncommon (e.g., 
iNice fc Tavlorl Il995t) . and implies that the millisecond 
pulsar period immediately after spin-up was close to its 
present day value. It seems likely that mil lisecond puls ars 
form with periods of a few milliseconds ( Backer 1998). 

The surface magnetic field strength of the pulsar 
according to conventional assumptions is Bo = 3.2 x 
10 19 (PP int ) 1/2 G = 1.9 x 10 8 G. This value is typical 
of millisecond pulsars. 

5.5. Timing Noise 

The PSR J1713±0747 timing data show timing noise. 
Analysis of the timing noise is complicated by the large 
gap in the data between 1994 and 1998 and because of 
the need to allow an offset in the arrival times before and 
after the gap. To analyze timing noise, we calculated 
residuals for a timing model fitting for only the pulsar 
spin-down parameters vq and v\ and an arbitrary offset 
between the pre-1994 and post-1998 data. Astrometric 
and binary parameters were held fixed at the values 
derived from the whitened timing model fit. The results 
are shown in Figure^,. 

Timing noise can be quantified by the fractional 
stability statistic a z , an adaptation of the "Allan 
variance" stat istic with modificati ons appropriate for 
pulsar timing ( Matsak is et al.lll9 97). In essence, a z for 
time interval r is calculated by dividing the post-fit 
residual arrival times, r(t), into intervals spanning length 
r, and fitting third order polynomials over each interval, 
r(t) — Y,3 =Q Ci(t — to) 1 . The first three terms of the 
polynomial are accounted for by the pulsar spin-down 
model. The final term, C3, is a measure of timing noise. 
The statistic is calculated by appropriately scaling the 
root-mean-square value of C3 from all intervals of a given 

length: a z (r) = (1/2^)^ (c 2 ) 1/2 . 

FigureElshows a z (r) of PSR J1713+0747 as calculated 
from the residual pulse arrival times shown in Figure^. 
For uniformly sampled data exhibiting white noise, o~ z 
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would fall off as r -3 / 2 ; this is the behavior observed on 
all but the longest time scale. At t — 12 yr, timing 
noise is evident. Because of the arbitrary offset fit 
between the 1994 and 1998 data, it is possible that the 
calculated value of o~ z (j) underestimates its true value 
on the longest time scale. 

The physi cal mechanism underlying timing noise is 
not known. lArzoumanian et al.l l)1994|) studied a large 
collection of pulsars using a statistic, A 8 , which is the 
logarith m of a scaled version of cr z (r) at r = 10 s s « 3 yr 
(see also lBackerll2004l) . From Figure [5J it is evident that 
at a time scale of 3yr, o~ z is dominated by measurement 
noise, so the measured values are upper limits on the 
intrinsic irregularities in the pulsar signal. Estimating 
logcr z (3yr) < -14.2 for PSR J171 3 +074 7, we calculate 
A 8 < —5.6. lArzoumanian et all Jj994) found A 8 to 
be correlated with pulse period derivative according to 
Ag = 6.6 + 0.6 log P. This formula predicts Ag = —5.5 
for PSR J1713+0747, close to the observed upper limit. 
Thus the rotational stability of J1713+0747 is as good 
(or better) than expected, despite the long-term timing 
noise. 

5.6. Solar System Ephemerides 

As discussed above, we used the DE405 solar system 
ephemeris to reduce the pulse arrival times to the 
solar system barycenter. To reduce pulse arrival times 
measured with uncertainty of only ~ 200 ns, the Earth's 
position with respect to the barycenter must be known 
with precision 60 m. We analyzed our data with both 
DE405 and its predecessor, DE200. To compare the 
two ephemerides, we performed timing fits on the the 
post-1998 data. Excluding the earlier data from these 
tests allowed us to avoid problems stemming from the 
arbitrary offset between the earlier and later data. Using 
each ephemeris, we fit the data to a full timing model, 



fitting for all the standard parameters but not allowing 
pulse frequency derivatives above v\ (i.e., no timing noise 
terms). The results are shown in Figure ITU1 The DE405 
ephemeris clearly gives a better fit. The differences in 
timing quality are easily explained by differences in the 
ephemerides themselves, particularly the incorporation 
of improved measurements of outer planet masses into 
DE405 (E. M. Standish, private communication). 

5.7. Testing Strong Field Gravity 

PSR J1713+0747 is one of several long-orbit, 
low-eccentricity binary pulsars which have been used 
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1996: Wcx 2000). The relevant observational signatures 
depend on the orientations of the binary systems under 
study, something not usually known; hence, probabilistic 
arguments have been used to constrain equivalence 
principle violations based on observations of several 
pulsars. In contrast, we have established the orientation 
(i and ft) to PSR J1713+0747, and we have measured its 
distance as well. This allows us to set the first absolute 
limits both on violation of the Strong Equivalence 
Principle (SEP) and on the magnitude of the strong-field 
equivalent of the Parametrized Post-Newtonian (PPN) 
parameter . 

5.7.1. The Strong Equivalence Principle 

If the SEP were violated, objects with different 
fractional mass contributions from self-gravitation would 
fall differently in an external gravitational field. This is 
quantified by the parameter 77, defined through m g /mi = 
1+riEg/mc 2 , where m g , mi, and E g are the gravitational 
mass, inertial mass, and gravitational self-energy of 
a body. Nonzero values o f t? would resul t in the 
polarization of binary orbits l|Nordt vedtl 1 1 968[) . Lunar 
laser ranging experiments set limits on the polarization 
of the Earth-Moon orbit in the gravitation al field of the 
Sun, constraining \r)\ to be less than 0.001 ijDickev et alJ 
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Fig. 9. — Noise statistic <r z as a function of time scale, r. 
Gaussian noise would result in cr z ~ r~ 3 / 2 , parallel to the dashed 
line. The upturn at r > 2 yr is indicative of timing noise with a red 
spectrum. The a z values at the longest time scales may be affected 
by the presence of a 4-year gap in the data and by the need for 
an arbitrary offset between pre- and post-gap data; however, we 
believe the upward trend at r > 2 yr to be qualitatively correct. 
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Fig. 10. — Residual arrival times from timing models for 
1998-2004 data calculated using (a) the DE200 solar system 
ephemeris and (b) the DE405 solar system ephemeris. 
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119941: 1 Willi 12001ft . In the case of a pulsar-white-dwarf 
binary, the orbit would be polarized in the direction of 
the gravitational pull of the Galaxy. The parameter 
to be constrained, A, is similar to T), but without the 
requirement of linear dependence on E g /mc 2 . It is 
defined for an individual body, a, by (m 9 /mj) a = 1 + A a ; 
dynamics of a binary orbit depend o n the difference A 
. \ — A2 between the two objects (see IDamour fc Schaferl 
1991, for full details). 

The (small) "forced" eccentricity ep of the orbit 
induced by the SEP violation may be written as: 



\ e F\ = 7T 



Ag ± c 2 



2 FG (m 1 +m 2 )(27r/P b ) 2 ' 



(8) 



where g± is the projection of the Galactic gravitational 
field onto the orbital plane and, in General Relativity, 
F = 1. This forced eccentricity may be comparable 
in magnitude to the "natural" eccentricity en(t) of the 
system, which by definition rotates at the rate of the 
advance of periastron and may, at any time, be oriented 
in such a way as to nearly cancel the forced eccentr icity. 
Defini ng the angle between these two vectors as 9, iWexl 
( 1997) writes the inequality: 



\e F \<e^(9), CiW = 



l/sin( 



1 

-l/sin0 



9 G [0, tt/2) 
9 e [tt/2, 3tt/2] . 
9 e (37r/2,27r) 
(9) 

where e is the observed eccentricity. 

The proj ection of g onto the or bital plane can be 
written as (Damour & Schafcr 1991): 



211/2 



|gj_| = |g| [1 — (cos i cos A + sin i sin A sm(cf> g — fi)) ] 

(10) 

where <p g is the angle of the projection of g onto the plane 
of the sky (effectively a celestial position angle, defined 
north through cast as for fi), and A is the angle between g 
and the line from the pulsar to the Earth. The value of |g| 
can be determined from model s of the Galactic potential 
(e.g.. lKuiiken fc Gilmorel 1989ft and the Galactic rotation 
curve; for a pulsar out of the Plane of the Galaxy A is 
given by 

d 2 + R\ + z 2 



cos A = 



-Rl 



2 d ^JR\ + z 2 



where R\ = (R 2 + (dcosb) 2 — 2 R a d cos b cos I) 1 ' ' 2 is the 
Galactic radius of the pulsar, R$ is the distance from the 
Earth to the Galactic center, z is the distance from the 
pulsar to the Galactic plane, and / and b are the Galactic 
coordinates of the pulsar. 

Historically both 9 and fl have been completely 
unknown for the pulsars used to test for violation of 
the SEP, and the tests have therefore made statistical 
arguments, assuming both angles to be uniformly 
distributed between and 2ir (e.g.. iDamour fc S chafcr 
IT991 . The masses have also been poorly constrained 
and thus av erages ove r likely model populations were 
also needed (Wcx 2000). Following these procedures, an 
ensemble of long-orbit pulsars yields a limit of |A| < 
0.009 at 95% confidence. 

With both masses and the angle of the line of 
nodes well-constrained for PSR J1713+0747, a single, 
more robust limit on violation of the SEP becomes 
possible. Because of the complicated dependence of 



A on the pulsar distance, as well as the asymmetric 
distribution of allowed masses, we obtain the limit via 
Monte Carlo simulation, assuming the parallax to be 
normally distributed and 9 to be uniformly distributed, 
and sampling the allowed m2-cos i-fl range according 
the probability distribution derived from the % S r id 
discussed above. We find that |A| < 0.013 at 95% 
confidence. This is nearly as good a limit as could 
be expected from this pulsar, as just over 90% of the 
Galactic acceleration vector is parallel to the plane of the 
orbit. The most stringent test of SEP violation, however, 
continues to rely on an ensemble of pulsars. 

5.7.2. Post-Newtonian Parameter 0:3 

The parameter is one of ten PPN parameters 
formulated to describe departure s from General 
Relat ivity in the weak-field limit l)Will fc Nordtvedil 
119721) . Its strong- field analog, da, can be tested 
by pulsar timing llDamour fc Esposito-Faresd Il992t 
IDamour fc Esposito-Fareselll996ft . A non-zero at would 
imply both violation of local Lorent z inyarianc e and 
non-conservation of momentum (e.g.. IWill 120011. For 
a pulsar-white-dwarf binary, the net effect would be 
to cause an accelerati on of the binary system given by 
jBfill fc namoiirilTflflffr 



a,-. 



. ; - gO^CpW x f2 Sp , 



(11) 



where fisp is the spin angular frequency of the 
pulsar (|f2s P | = 2ttuq), w is the absolute velocity 
of the system, and c p denotes the compactness 
of the pulsar, roughly the fraction of its mass 
contained in gravitational self-energy. An approximate 



expression for the compactness i s 



0.21toi/M q 



ijDamour fc Esposito-Faresel 11992ft . As above, this 
acceleration will induce a "forced" eccentricity in the 
orbit, given by 



a 3 



=p|w| P b 2 



24tt P G(mi +m 2 ) 



sin/3, 



(12) 



where (3 is the angle between w and f2s p , an d P 1S 
the spin period of the pulsar: P = 27r/f2s p . Using an 
ensemble of pulsars and as suming rand om distributions 
of the orbital inclinations, IWexl <|2000f1 found a limit of 
|&3| < 1.5 x 10~ 19 . To calculate the limit imposed by 
PSR J1713+0747, we again proceed via Monte Carlo 
simulation, sampling the parameters as described above. 
We assume that the pulsar spin and orbital angular 
momenta have been aligned during the spin-up episode 
that recycled the pulsar. Thus with the full orientation of 
the orbit known for each simulation point, we can easily 
calculate (3 given any w. This absolute velocity is taken 
to be in the reference frame of the Cosmic Microwave 
Background (CMB), and is calculated from the motion of 
the Solar System in the CMB reference frame (369 km/s 
toward a = 11^3, 5 = -7°; iFixsen et aDl!994j1 and the 
three dimensional velocity of the binary system relative 
to the Solar System. The radial component of this 
second velocity vector is unknown. We checked a range 
of radial velocities from —200 km/s and +200 km/s; for 
each simulation point, we adopted the value that gives 
the smallest projection of the total velocity onto the 
plane of the orbit. This value was always in the range 



11 



—60 to +25km/s. We thus arrived at an absolute 
95%-confidence limit of |d 3 | < 1.2 x 1CT 19 . This is better 
than the limit derived from the ensemble of pulsars and, 
as it is less statistical in nature, it may be considered a 
robust true limit on d^. Future refinement of the timing 
parameters will improve this limit somewhat, with the 
floor ultimately determined by the orbital geometry. 

6. CONCLUSION 

PSR J1713+0747 has lived up to its promise as one 
of the best pulsars for high precision timing, with 200 ns 
residual pulse arrival times attained on time scales of 
years, and residuals well under 2/is over the full 12-yr 
data set. Major findings include: 

1. The pulsar mass is constrained to mi = 1.3±0.2 M Q 
by the measured Shapiro delay. If the secondary mass is 
restricted to values predicted by the theoretical orbital 
period-core mass relation, the pulsar mass is somewhat 
higher, To! = 1.53T°;°*. 

2. The parallax is 7r = 0.89 ± 0.08 mas, corresponding 
to a distance of 1.1 ± 0.1 kpc. This is consistent with 
predictions based on the pulsar's dispersion measure. 

3. The orientation of the binary has been fully 
determined by the combined measurements of Shapiro 
delay and of perturbations of orbital elements due to 
relative Earth-pulsar motion. The orientation and very 
low ellipticity of the orbit lead to an improved constraint 



on deviations from general relativity. 

The timing precision attainable on short time scales 
appears to be limited by measurement precision. There 
is every reason to expect timing precision to improve 
in the coming years as a new generation of wide 
bandwidth coherent dedispersion systems are employed 
at radio telescopes, directly resulting in more precise 
measurement of the Shapiro delay and hence of the pulsar 
and white dwarf masses. 
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